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On example of diffusion-limited reversible A + A ^ B reactions we re-examine 
two fundamental concepts of classical chemical kinetics - the notion of "Chemical 
Equilibrium" and the "Law of Mass Action". We consider a general model with 
distance-dependent reaction rates, such that any pair of A particles, performing 
standard random walks on sites of a d-dimensional lattice and being at a distance 
(i apart of each other at time moment t, may associate forming a B particle at 
the rate k + {^). In turn, any randomly moving B particle may spontaneously dis- 
sociate at the rate k-{\) into a geminate pair of As "born" at a distance A apart 
of each other. Within a formally exact approach based on Gardiner's Poisson rep- 
resentation method we show that the asymptotic t = oo state attained by such 
diffusion-limited reactions is generally not a true thermodynamic equilibrium, but 
rather a non-equilibrium steady-state, and that the Law of Mass Action is invalid. 
The classical concepts hold only in case when the ratio k + (fi)/k-(fi) does not depend 
on fi for any fi. 

PACS numbers: 05.70.Ln; 05.40.-a; 05.45.-a; 82.20.-w 
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I. INTRODUCTION 

"Chemical Equilibrium" (CE) and the "Law of Mass Action" (LMA) are two central 
concepts of classical chemical kinetics (see, e.g., Refs. In virtually every text-book 

one finds, regarding, for instance, the behavior of reversible association/dissociation reaction 
of the form 

A + A S B, (1) 

K- 

where K + and K_ are the forward and the backward rate constants, respectively, that the 
asymptotic state achieved in closed systems at t = oo is the state of Chemical Equilibrium - 
state with no net change in activity, or concentration with time t. Thermodynamically, the 
condition of CE is the condition in which the driving forces of the reaction in Eq. (0) (or any 
other reversible reaction) are equal and opposite. This condition implies that no spontaneous 
change is observed and that, according to the zeroth principle of thermodynamics, the net 
Gibbs free energy change of a mixture of reactants and products vanishes. Kinetically, 
the condition of CE is the condition in which the rates of the forward, K+a^, and the 
backward, K-boo, reactions are equal and opposite, such that Oqq and b^ - the "equilibrium" 
concentrations of A and B species, obey the LMA: 

-K+al + = or ^ = = K eq , (2) 

with K eq being the "equilibrium" constant, dependent only on the thermodynamic properties 
of the reactive system . It is important to remark that once the CE is achieved, the 
forward and backward reactions continue to run. It is just at equilibrium, since the rates 
are equal, there is no visible or measurable change in the system. 

In this paper we re-examine these two fundamental concepts, serving to define the compo- 
sition of reactive mixtures (as well as general trends in reactive systems) , in case of reversible 



reactions w 
reactions 
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ve diffusive particles and products, i.e. the so-called diffusion-limited 



Our focal questions here are whether for such diffusion-limited 
reactions taking place in closed systems the CE is always a true thermodynamic equilib- 
rium state (TES), and whether the LMA in Eq.(|2*|) always holds. We emphasize that it 
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does not mean that the principles of thermodynamics are contested: according to the zeroth 
principle, for a closed system the CE always exists and is by definition the thermodynamic 
equilibrium state. What is meant here by true thermodynamic equilibrium state, is much 
more restrictive: within the conventional picture it designates a state fully described by 
the thermodynamical, as opposed to dynamical, quantities, and which satisfies the detailed 
balance equilibrium. The fact that the CE is not, in some cases described below, a TES 
is therefore fully compatible with the laws of thermodynamics. Note that in particular the 
LMA is not a law of thermodynamics, since it relies on an ideal gas approximation which is 
not necessarily always the case. 

We concentrate on a particular reaction scheme - the simple association/dissociation re- 
action in Eq. ([T|) . but our analysis can be readily generalized for any other type of reversible 
reaction. We consider here a rather general lattice model of reactions in Eq.(^), which was 
first studied analytically by Zeldovich and Ovchinnikov [18]. In this model the A particles 
and products B perform standard random walks on sites of a d- dimensional hypercubic lat- 
tice and the elementary reaction rates are long-ranged and dependent on the instantaneous 
distance between any two A particles; that is, any pair of As may associate (forming a 
B particle) at any moment of time t at rate k + (fi), where fi is the instantaneous distance 
separating these two particles. In turn, any B particle may spontaneously dissociate at 
rate k_(X) giving birth to a geminate pair of As separated by a distance A. In our analy- 
sis we suppose that the bimolecular elementary reaction rate k + (fi) and the unimolecular 
elementary reaction rate fc-(A) are arbitrary (integrable) functions of [i and A. Note also 
that we have chosen the lattice formulation just for the convenience of exposition; an analo- 
gous continuous-space formulation can be readily worked out. For this model, in terms of a 
formally exact approach based on Gardiner's Poisson representation method P|, we obtain 
exact non-linear Langevin equations describing the time evolution of complex-valued Poisson 
fields, whose mean values determine the A and B particles' mean concentrations. Solutions 
of these Langevin equations in the asymptotic t = oo state are obtained by two different 
approaches: via a) a certain decoupling approximation and b) a systematic diagrammatic 
expansion. From these solutions, which coincide in the leading order, we deduce a general 
criterion determining the conditions when the classical LMA in Eq.(0) holds and when the 
asymptotic t = oo state is a true thermodynamic equilibrium. We show that this may only 
happen when the distance-dependent elementary reaction rates obey a rather strong (and 
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apparently unrealistic) condition: the ratio k + (/i) / '&_(//) does not depend on /i for any //! In 
case when this microscopic restriction is violated even at a single point, one can show that 
the detailed balance is broken. Here we demonstrate that this violation of the detailed bal- 
ance at the microscopic scale has macroscopic consequences: the LMA in Eq.(J2J) is violated, 
particles' concentrations are spatially correlated, the correlation length is macroscopically 
large, and, remarkably, the CE is not a true thermodynamic equilibrium but rather a non- 
equilibrium steady-state, depending on dynamical properties such as the particles' diffusion 
coefficient. The profound reason of these spectacular macroscopic effects is that the break- 
ing of detailed balance generates a non-vanishing probability current which modifies and 
sustains the fluctuations. In turn these fluctuations correlate the particles' concentrations, 
perturbate the rate of the forward reaction and thus displace the equilibrium concentrations, 
breaking the LMA. Note that such a current occurs in the phase space and there is no net 
transport in the real space. We remark that this non-equilibrium steady-state breaking the 
detailed balance provides an example of irreversible circulation of fluctuations - a notion 
put forward in Ref.[10]. Using system size expansion method, Tomita and Tomita have 
demonstrated that for statistical physics systems a non-vanishing probability current can be 
generated by breaking of the detailed balance (a state which they call a "cyclic balance"). 
We note parenthetically that appearance of a non-equilibrium steady state is generic for 



lll | and subsequently elab- 



open reaction-diffusion systems, as first suggested by Kuramoto 

orated in Refs. 3,0, 0,0]- Under certain conditions, it may also take place for reactions 
in which the particles number is not explicitly conserved [16]. We emphasize that here such 
a non-equilibrium steady-state emerges in a closed system with strictly conserved overall 
concentration of particles and products, without any external inflow of particles! This gives 
a striking example of a steady state breaking the detailed balance equilibrium, characterized 
by non vanishing probability currents in the phase space. 

We note that despite the common belief that in closed systems the CE is a true TES and 
that the LMA may be taken for granted, there are some good reasons which might make 
our questions legitimate. Indeed, the classical kinetic picture has already been proven to be 
inadequate in many situations and a number of significant deviations from the text-book 
behavior has been discovered. These deviations concern primarily the kinetic behavior. 
For example, for reversible reactions conventional chemical kinetics predicts an exponential 
approach toward the CE state. It has been realized, however, that this is not the case 
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for reversible reactions involving diffusive species; here, the concentrations approach the 
asymptotic state t — > oo only as a power law {t~ d l 2 in d dimensions)! Such a behavior stems 
out of many particle and non-linear effects. Because of the bimolecular (non-linear) forward 
reaction, the time evolution of observables - particles mean concentrations, appears to be 
coupled to the evolution of the pairwise correlation functions. In turn, long-time decay of 
correlations is dominated by existence of the conserved values, e.g., the overall concentration 
of particles and products, which are not affected by reactions and thus represent purely 
diffusive modes of the system. Theoretically, the power law approach to the asymptotic 
state has been first predicted in Ref . \v\ using physical arguments based on the analysis of 
spatial concentration fluctuations, and subsequently elucidated in terms of more elaborated 



approaches 
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solutions have been obtained 
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] . In some cases, exact dynamical many-particle 
As well, this power law behavior has been indeed 
observed in excited state proton transfer reactions 3(], 31 1, and also seen in numerical MC 
simulations \?>\ - 

The question whether in closed systems the LMA in Eq.(j2J) is valid for the reversible 
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diffusion-limited reactions at t = oo has already been raised in Refs.|18| and 
Refs. ^| an approximate approach has been proposed to describe kinetics of reversible asso- 
ciation/dissociation reactions in Eq.fJJ) with distance-dependent elementary reaction rates. 
In this approach, the evolution of the reactive system has been reduced to the evolution of a 
mixture of two quantum Bose gases, which was solved by extraction of the condensate and 
approximate second quantization method, similar to the Bogolyubov's theory of a weakly 
non-ideal Bose gas In Refs. , some corrections to the LMA non- vanishing in the limit 
t — > oo have been found. Corrections to the LMA have also been obtained in Ref.Q in 
terms of a suitably extended Smoluchowski approach for reversible, microscopically inhomo- 
geneous contact reactions, such that the reaction radius R of the forward reaction is unequal 
to the radius A of pairs appearing within the course of the backward reaction. In three di- 
mensions, it was found that in the asymptotic t = oo state the particles concentration obey 



where D is particles diffusion coefficient. Note that according to Eq.(|3*j). particles concen- 
trations in the asymptotic t = oo state are dependent on such "kinetic" parameter as D, 
which could be quite alarming, if Eq.Q were derived in terms of a more reliable approach. 
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Further on, Refs.|22j have pursued a different uncontrollable approach analyzing the tem- 
poral evolution of several contact diffusion-limited reactions in terms of reaction-diffusion 
equations for local concentrations. In these equations, the stochastic nature of particles 
transport and reactions have been incorporated by random source terms, derived within the 
framework of a hydrodynamic-level stochastic description of number densities for a dilute, 
chemically reacting system Ref . [34| . Solving the resulting non- linear Langevin equations un- 
der an assumption that the concentration fields are Gaussian (which implies an automatic 
decoupling of fourth-order correlations and hence, insures that the third-order correlations 
vanish), Refs. 22j predicted not only non- vanishing corrections to the LMA, but also shown 
that in the asymptotic t = oo state particles' spatial distributions are correlated. Moreover, 
it has been realized that the correlation length depends on particles' diffusion coefficients, 
which rules out, of course, that the CE is a true TES, but rather represents a non-equilibrium 
steady- state. 

On the other hmd , formal* exact M on S ohtahred in Refsfl aad Q fa severa, 
reversible diffusion-limited reactions have demonstrated that the LMA in Eq.(J2J) is strictly 
valid at t = oo, and that the asymptotic t = oo state has a Poissonian spatial distribution 
of concentrations, which sig nifies that the CE is a true TES. Consequently, approximate 
^> 21, 3 an d rigorous |2tL approaches, although agree on the dynamical behavior, are 
at some odds concerning the properties of the asymptotic t = oo state. 

Both exact approaches (23, l2i| > however, focused on contact diffusion-limited reactions 
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with some rather restrictive conditions imposed on the elementary reaction acts. More specif- 
ically, in Ref. (22]], which considered a continuous-space model of different reversible reactions 
and made use of suitably generalized quantum field theory techniques of Refs. (35, 
forward reaction radius and the radius of the geminate pair born in the backward reaction 
act were both set equal to zero. In Ref.j^J, which used a lattice formulation ingeniously 
extending the Poisson representation method of Gardiner [9], it was stipulated that the for- 
ward reaction takes place when both species occupy simultaneously the same lattice site, 
while the dissociation of the reaction product produces a geminate pair of reactants born 
at the same lattice site. We note parenthetically that for this particular case both Refs. 
and [2l|, Eq. (|HJ) , predict that corrections to the LMA vanish as t — > 00 and thus the LMA 
obeys its classical form in Eq.(j2*|). 

Strictly speaking, reaction events are not describable within the framework of a classical 
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theory only; an elementary reaction act results from an interplay of many factors and is 
influenced by solvent structure, potential interactions, a variety of particles' ener gies and 
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angular orientations, quantum processes of different origin and etc [ 
At such scales, a notion of a fixed "reaction radius" R does not make much sense - a 
"reactive" boundary condition, imposed at a fixed distance R separating the reactants, 
is just a mathematical trick employed to obtain a tractable formalism and to circumvent 
enormous technical difficulties. This technical breakthrough is achieved, however, at expense 
of introducing a certain degree of arbitrariness regarding the choice of the value of R and 
of the reaction rate itself, which both, in consequence, are rather ill-defined. In this regard, 
for a more adequate description of an elementary reaction act, (still being, however, in a 
reasonable compromise between either too restrictive or too complicated theoretical model) 
it is appropriate to introduce, as it has been done, in particular, in Ref. 3|, a distance \x 
dependent elementary reaction rates for both forward and backward reactions. Note that 
an account for the long-range distance-dependent character of elementary reaction rates is 
indispensable for the analysis of kinetics of reaction including remote electron or proton 



transfer |4j, |4Jjj , for which the contact approximation is meaningless. It appears, as well, 
that it is indispensable also in the general case for the analysis of such a delicate issue as the 
nature of the asymptotic t = oo state. The meaning of this statement will become clearer 
as we proceed. 

This paper is outlined as follows. In section 2 we formulate the model, write down the 
master equation describing particles reactions and diffusion on the microscopic many parti- 
cle level, and, within the framework of Gardiner's Poisson representation method jjj, derive 
non-linear Langevin equations describing the time evolution of the Poisson fields. In sec- 
tion 3 we discuss solutions of the non-linear Langevin equations using a certain decoupling 
approximation. Within this approach, we define a criterion determining the conditions 
when the classical LMA in Eq.(j2| holds and when the asymptotic t = oo state is a true 
thermodynamic equilibrium. Focusing next on a particular case of microscopically inhomo- 
geneous contact reactions, we determine explicitly corrections to the LMA. Further on, in 
section 4, we set up a systematic approximation scheme previously developed for contact 
reactions in Ref. [27J • We show, within this mathematically rigorous approach, that the cri- 
terion obtained in section 3 is exact. Moreover, we demonstrate that for microscopically 
inhomogeneous contact reactions, the corrections to the LMA obtained in section 3 within 
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the decoupling approximation, as well as the result in Eq.(JH} obtained in Ref. 2l( within a 
suitably extended Smoluchowski approach, are exact in the linear order in deviation from 
the equilibrium. Besides, we present here several explicit results for exponential reaction 
rates. Finally, we conclude in section 5 with a brief summary of results and discussion. 



II. MODEL AND EVOLUTION EQUATIONS 

Consider an infinite <i-dimensional hypercubic lattice of spacing £, containing particles of 
two types - A and B, which perform unconstrained (an arbitrary number of particles can 
occupy any lattice site) random walks between neighboring sites. Any two A particles, being 
at a distance fi from each other, may enter into reaction at the rate k + (n) forming a single 
B particle, placed at the half-distance between two As. Further on, any B particle at any 
moment of time may spontaneously dissociate at the rate k-(X) producing a geminate pair 
of A particles born (with a random orientation) at a distance A apart of each other. Note 
that the true reaction constants, (those entering Eq.(J2J)) are determined as K + = k + (fi) 
and K_ = ^ A fc_(A) Most of our analysis will be performed supposing that k + (fi) and 
fc-(A) are arbitrary integrable functions of fi and A. 

In what follows, we will distinguish between two situations: 

(i) the case of microscopically homogeneous elementary reactions, when k + (fi) is strictly 
proportional to k-(fi) (or, in other words, when the ratio k + (fi) / k-(fi) is independent of \i 
for any /i) 

(ii) the general case of microscopically inhomogeneous reactions, where k + (/i) / k_(fi) is /i- 
dependent, at least in some region of space. 

The criterium of microscopically homogeneous reactions has an important interpretation. 
Indeed, for microscopically inhomogeneous reactions, any steady distribution of the concen- 
trations of the species A and B breaks the detailed balance equilibrium. This is due to the 
fact that in this case the equilibrium constant /c_(/z)/fc + (/z) is space dependent : equilibrium 
can not be satisfied at each point for homogeneous concentrations. We will show below that 
this violation of the detailed balance has important implications at the macroscopic scale. 

The state of the system at time t is determined by the time-dependent numbers A(x) and 
B(x) of A and B particles at site x of the lattice. The set of such numbers is denoted as {A} 
and {B} and P[A,B,t] stands for the probability of finding the system at time moment t 



9 



in the {A} and {B} state. 

Our analytical approach is based on the formally exact Poisson representation method, 
proposed originally by Gardiner |9'] , and subsequently generalized in Ref . J3] for the descrip- 
tion of the fluctuation-induced kinetics of reversible, microscopically homogeneous contact 
diffusion-limited reactions. Extension of this approach to the case of distance-dependent 
reactions is straightforward and here we merely outline the steps involved. 

The starting point of our analysis is the following master equation describing the time 
evolution of P[A, B, t]: 

d t P[A 1 B;t] = CP[A,B;t] } (4) 

where £ is an operator, accounting for the reaction and diffusion processes on the micro- 
scopic, many particle level. Explicitly, the right-hand-side of Eq. (JH) of reads: 



D 



D 



CP[A,B;t] = ™^E \[A(e x ) + l)P[A(x)-l,A(e x ) + l,B,t}-A(x)P[A,B,t] 



+ p E E [{ B ^ + l ) P ^ B & - B ^ + M - B(x)P[A, B, t] 

x e x 

+ E k ^ E [{ A ( x - f ) + 1 + <w) ( A ( x + f ) + 1 



x P[A(x - |) + 1 + 6^ , A(x + ^) + l + 6^ , B(x) - 1, t] 
- A(x-^){A{x + ^)-8^)P[A,B,t] 

+ J2 k -^) E ( B ( x ) + l ) P i A ( x - f ) - 1 - <W>, A ( x + f ) - 1 - <W, + 1, t] 

fl x 

where D = £ 2 /2dr is A and 5 particles diffusion coefficient, r is a characteristic hopping 
time, while the symbol denotes summation over all possible orientations of the lattice 
vector e x . As the initial condition to Eq.©, we choose an uncorrelated Poisson distribution 
on each site and for each species: 



A A(x) R B(x) 

P[A, B,t = 0} = e- A °- B « TT 2 "°. - , 
L 1 11 A(x)\B(x)\ 

where Aq and Bq are mean occupation numbers of A and B particles per lattice site. 



(6) 
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The next step consists in projecting the particles numbers on sites x onto the Poisson 
states jjj: 



P[A,B,t) = J Y[da(x)dp( 



, exp (-a(x)) a(x) A( - x ^ 

c) ip x 



x exp (-P(x))P( x )" { " ] F[a(x),p(x),t] (7) 

JD[X). p. 

Existence of such a transformation (and of the inverse one) is well established |9|. As a 
matter of fact, the factorial moments of particles' numbers, say of A(x), are related to the 
Poisson fields a(x) and j3(x) via a simple formula: 

(A(x)) n = (A(x)(A(x) - l)...(A(x) - (n- 1))) = 
= / Y[da(x)d(3(x) a n (x)F[a(x),p(x),t] = (a n (x)) , (8) 

where the angle brackets denote the probabilistic average stemming out of stochastic reaction 
processes. Note, however, that the Poisson fields a(x) and (3{x) may be complex- valued; 
consequently, F[a(x), 0(x), t] may admit negative values and hence can not be interpreted 
CIS db probability distribution. 

Defining next the generating function 

G[u{x),v(x),t] = J2H u{x) a v{x) b P[A, B, t] (9) 

a,b x 

one represents G[u(x),v(x),t] as 

G[u(x), v(x), t] = / Y\ da{x)df3{x) exp[(w(x) — l)a(x)] x 

X 

x exp[(v(x) - l)/3(x)]F[a{x),l3(x),i\, (10) 

which yields, eventually, the following Fokker-Planck equation determining evolution of 

F[a(x),(3(x),t]: 

r 9 _ r nl d 



d t F[a(x),f3(x),t] = -EEl^J^^^ + ^y^/ 3 ]- 

-C 3 [a,(3])F[a(x),(3(x),t] (11) 



X /i 

_ i d?_ 

2 da{x - n/2)da{x + /i/2) 

where 

C±[a, P] = —k + (fi)a(x) (a(x — /jl) + a(x + + k-{fi) (j3(x — fj,) + (3{x + /i) 
< C 2 [a,(3] = k+(n)a(x - fj,/2)a(x + fi/2) -k-(n)p(x) + DAp, 
C 3 [a, (3] = k + ([i)a(x + fi/2)a(x — /x/2) — k_(p,)(3(x), 

(12) 



+ DAa, 
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where A is the second finite difference operator. 

Note that Eq. fjllj) can not be, of course, solved exactly. On the other hand, we are in 
position to draw one important conclusion directly from the form of Eq. ffllj) . By inspection, 
one may verify that because of the presence of "mixed" derivatives d 2 /da(x — fi/2)da(x + 
fi,/2)..., in the microscopically inhomogeneous case, when k+{p) / k-(fi) is /i-dependent, this 
Fokker-Planck equation does not admit any stationary solution of the form: 

F[a(x),P(x),t] = J] 5[a(x) - (a)]5[(3(x) - ((3)} (13) 

X 

Therefore, if k + (^)/k-(fi) depends on //, the stationary distribution P[A, B, t] is not a prod- 
uct of uncorrelated Poisson distributions, as one obtains in case of contact, microscopically 
homogeneous reactions. This implies, in turn, that the detailed balanced is not satisfied, 
which hints us that here the asymptotic t — oo state might be rather particular. Indeed, be- 
low we will demonstrate that the fact that F[a(x), /3(x), t] does not factorize into the product 
EL fl a ( x )i P( x )i t]i where each multiplier f[a(x), /3(x), t] depends only on the Poisson fields 
a(x) and f3(x) on the site x, implies that the CE state is not a true equilibrium state and 
that the pairwise correlations in particles' spatial distribution exist, which, in consequence, 
results in violation of the LMA in Eq.Q. 

Using Ito's equivalence, we find the following Langevin equations corresponding to equa- 
tion (fTTl): 

d t a(x) = ^2Ci[a,0\ + COM), 

A ( 14 ) 

d t fi(x) = ^2C a [a,0\, 

where ({x,t) is a Gaussian noise with zero mean value, whose correlation is given by 

(C(x,t)C(x + x',t')) = 5(t-t')[k4x')((3(x')) 

- k+{x')(a{x-x'/2)a{x + x'/2)^ . (15) 

Before we proceed further, some comments are in order. We first remark that exactly 
the same Langevin equations can be obtained extending to the general case with distance 
dependent elementary reaction rates the field theoretical technique, developed in Ref. 27 
for contact microscopically homogeneous reactions. We believe, however, that the approach 
used here is more simple and transparent. Second, we note that contrary to the situation 
with contact microscopically homogeneous reactions (see Eqs.(12) and (13) in Ref. m 
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the general case when the ratio k + (n) / k_(n) is //-dependent, the amplitude of the noise 
correlation can not be expressed as a time derivative of < a(x) >; hence, noise correlation 
does not vanish as t — > oo. Finally, we note that non-linear Langevin equations have been 
already used in Refs. I21I I22I ] to analyze the temporal evolution of several contact diffusion- 
limited reactions. Noise terms, employed in this approach, were previously derived within the 
hydrodynamic- level stochastic approach in Ref. 34j . Our Eqs. (jl4|) and ()15|) have essentially 



the same structure as those used in Refs. [21! |22| with two notable exceptions: 
(i) Eqs.flHJ) and (JT3J) are obeyed by the Poisson fields, which are not physical concentrations, 
and moreover, the noise is complex-valued, while the corresponding Langevin equations in 
Refs. 121, 22] are formulated for the true local concentrations. 



(ii) the noise terms used in Refs. 



21 



221 ] contain an additional, compared to our Eq. (|15jl . 



term ~ DV ' x A(x)V ' x 8{%) , which does not vanish as t — > 00. This is apparently incorrect and 
represents an evident shortcoming of the hydrodynamic-level stochastic approach in Ref.Q]. 
Moreover, this is precisely the reason why the deviations from the LMA were obtained in 
Refs. |2l|, I22J for microscopically homogeneous contact reactions. 



III. THE ASYMPTOTIC t = 00 STATE: DECOUPLING APPROXIMATION. 

For convenience, we turn from now on to the continuous-space limit (£ — > 0), consider 
x as a continuous variable and operate with particles' concentrations a(x) = A(x)/l d and 
b(x) = B(x)/l d , instead of particles' numbers. Note also that for an accurate "passage" to 
such a limit, one has to introduce the rescaled reaction rates k + — > k + £ d ~ 1 and — > 
noise ( — > (/l d , as well as the Poisson fields a(x) — > a(x)/l d and (3(x) — > j3(x)/l d j^ . 

Taking into account that, in virtue of Eq.(jHJ, we have < a(x) >=< a(x) >= a and 
< (3{x) >=< b(x) >= b, we write the Poisson fields as follows: 

a(x, t) = a + Sa(x, t), 
/3(x,t) = b + 5/3(x,t) 

where a and b are (time-dependent) particles' mean concentrations, while 5a(x,t) and 
5/3(x,t) denote local deviations of the Poisson fields from their mean values. Hence, by 
definition, < 8a(x, t) >=< 5(3(x,t) >= 0. 

Substituting next expressions in Eq. ([TK|) into Langevin Eqs.((T3J), averaging the resulting 
equations and using the definition of the overall reaction constants K± = f k±(fi) dfi, we 
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find that particles' mean concentrations obey: 



d t a 



d t b 



-2K + a 2 + 2K_b + 2Q(0), 
K + a 2 - K-b-Q(0), 



(17) 



where 



fi(0) = - / dfik + ([i)a aa (fi), 
and a aa (/i) is the pairwise correlation function, 

^(/i) = (5a(x — n/2, t)5a(x + ji/2, t) 



(18) 



(19) 



of the fluctuations in the Poisson fields. The corresponding correlation function for the fluc- 
tuations in particles concentrations in the asymptotic t = oo state can be readily expressed 
in terms of cr aa (/i): 



c(/i)oo = ( ctoo - a{x, t = oo) ) ( cioo - a(x + /i, t = oo) ) ) = a^n)^ + a^S^) (20) 



Note that Eqs. (jl7|) are formally exact for any t. These equations show that the time 
evolution of observables - the particles' mean concentrations, is ostensibly coupled to the 
evolution of pairwise correlations in the reaction-diffusion system under study. 

Turning next to the infinite time limit, t — > oo, we get, slightly rearranging Eqs. (jl7j) . that 
the particles' mean concentrations obey: 



where the subscript oo signifies that we deal with the asymptotic t = oo solution. 

Note that Eq.([2ip resembles the classical Law of Mass Action in Eq.(J2J), but differs from 
it due to an extra term fi(0)/i^ + 6oo, which is dependent on the pairwise correlation function 
0"aa(/-Ooo and embodies all non-trivial physics associated with the fluctuation effects. There- 
fore, one may judge directly from Eq.(|21jl that the classical LMA holds only if A particles 
distribution is spatially uncorrelated in the asymptotic t = oo state, i.e. o" aQ (/i)oo = 0. 

Now, our aim is to determine the pair correlation function a aa (fi). This turns out to 
be, however, quite a complicated problem. More specifically, when one writes the evolution 
equations obeyed by pairwise correlations, he gets that these are coupled to the third-order 
correlations and etc, which is a signature of a genuine many particle problem. 







(21) 
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In this section we resort to an approximate approach, while some exact results will be 
presented in the next section. Our approach here is based on the assumption that the 
deviations 5a(x,t) and 8(3(x,t) are Gaussian random fields, which implies an automatic 
decoupling of fourth-order correlations into the product of pairwise ones, and hence, trun- 
cation of the hierarchy of coupled reaction-diffusion equations on the level of third-order 



correlations. Such an approach has been first employed in Ref. 



to obtain the t d / 4 -law 



describing the fluctuation-induced kinetics of irreversible A + B — > reactions, and subse- 



quently generalized in Refs.|21| and 



22 



for other irreversible and reversible reactions. 
Within this approach, we find that pairwise correlations obey the following system of 
equations: 

d t cr aa (n) = 2DA fl a aa (fi) - 2a J dX k+(\) (2a aa (n) + a aa (fx - A) + a aa {n + A) ) 



+ 2 / dX fc_(A) [cTofliU - A) + (Tapfji + A) J + 

+ k-(fjL)b - k + (n)a 2 - k + (n)a aa (n), (22) 



dtPap(ji) = 2DA fl a a(3 ([i) - (2aK + + K^cr a/3 (fi) - a J dX k+(X) (a a p(n + A) + a a(3 (fi - A) 

/A A 
dX k+(\) (cr aa {[i + - ) + a aa (fjL - - jj + 

+ J dXk-{X)(a l30 {fJ l + X)+a m {fj l -X)^ (23) 

and 

/A A 
dX k+(X) (a a p(n + - ) + <Tap(n - -)) , (24) 

where the "a/3" and " (3(3" pair correlation functions are defined as 



o*p{p) = [Sa{x-fi/2,t)5(3{x + fi/2,t)j, (25) 

and 

appifj) = (6P(x-n/2,t)6/3(x + n/2,t)y (26) 
Further on, introducing a pair of Fourier transforms: 

/(p) = JdfiJMffa), and f^) = -±-j 'dpe-^fip), (27) 
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we get from Eqs.(|22J) to (jUJ) that a k i(p), (k,l = a, /?), obey: 

d t cr aa (p) = -2^Dp 2 + 2aK + + 2ak + (p) s ja aa (p) 

+ 4k4p)a aP (p) + k_(p)b - k + (p)a 2 + fi(p), (28) 
d t a a p(p) = -(2Dp 2 + 2aK + + AT_ + 2ak + {p))a a p{p) 

+ 2ak + (p/2)a aa (p) + 2fc_(p)cr^(p), (29) 
<V^(p) = 4afc+(p/2)0- o/J (p) - 2(/V + K^a^ip), (30) 

where we have assumed that fc+(/i) and depend only on the absolute value of /x, and 

denoted 

fi(p) - -J dfxe^k + (fx)a a M- (31) 

Turning now to the asymptotic limit t = oo, we get from Eqs. (j28j) to ()3()j) that 0" aQ (p)oo is 
determined as the solution of the following equation: 



where 



aa (p)oo = F{p) fi(p)oo - e(p)oo , (32) 

F( , = X 2 (p) (X^p) + X 2 (p)) - X 3 (p) 
{P) 2(X l (p)X 2 (p)-X 3 (p))(X l (p)+X 2 (p)y { ) 

Xi(p) = Dp 2 + 2K +aoo + 2k+(p)a 00 , 
X 2 (p) = Dp 2 + K„, 

X 3 (p) = 4fc+(p/2)A;_(p)a 0O , (34) 



and 



t{p) OQ = k + {p)a 2 00 -k_{p)b 00 . (35) 

Note that Eq.(|52*|) only implicitly defines u a a(p)oo, since the term Q(p)oo on the right-hand- 
side of this equation is itself dependent on 0" Q , Q (p) oo , Eq.([3"T|). Therefore, Eq. (j3^j) is an integral 
equation and its solution may be found only if we specify and k_(/i). 



A. Criterion of validity of the LMA and of the True Thermodynamic Equilibrium: 

Decoupling Approximation. 

Despite the fact that (T aa {p) oo has been so far only implicitly defined, and in order to get 
some explicit results we have to fix k+(p) and we are in position now to deduce a 
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general criterion specifying when the LMA in Eq.(J2J) holds and when it is violated. To do 
this, we have merely to define conditions when a aa (p)oo 7^ 0. 

We proceed as follows. First, we set in Eq.(jH2} o~aa(p)oo = (i.e. cr a a(/-0oo = 0). This 
implies that fi(p)oo = 0, and hence, in order to have a aQ (p)oo = 0, the following relation 
should hold as an identity: 

k-(p)boo = /c+(p) a L for any p, (36) 

or, in the /i-domain, 

/c_(/x)feoo = k + (p,)a 2 x> for any /i. (37) 

Note that in this only case the amplitude of the noise- noise correlation in Eq. (jl5j) is propor- 
tional to dta and hence, vanishes as t — > 00. 

Now, one readily verifies that the identities in Eqs. (j36|) and ()37|) may hold only if the 
reactions rates are identic functions of /1, i.e. the elementary reactions are microscopically 
homogeneous. If this condition is violated, even at least at a single point, pair correlations 
(Tea (/x)oo are non-zero, and hence, in virtue of Eq. (j21j) . the LMA in Eq.(J2J) is violated. 
Hence, we find that the conventional LMA in Eq.(|2"|) is always violated for microscopically 
inhomogeneous reactions. 

Next, for microscopically inhomogeneous reactions, we estimate the decay of pairwise 
correlations in the limit /x — > 00. In the limit p — > 0, (p — > 00), the kernel F(p) in Eq. (j33|) 
attains the form: 

2D V + Dp 2 (AK+aoc + 3K.) + K % _ 

F(p) ~ (38) 

2Dp 2 (Dp 2 + AK +aoo + K_j (2Dp 2 + AK +aoo + K_J 

Now, let A r be the characteristic decay length of k + (p). Since this property is usually of a 
range of a few interparticle separations, one may expect that cx aQ (/i) oo varies with /1 much 
more slowly than k + (fi). This assumption will be checked for consistency afterwards. Hence, 
we may approximate f2(p)oo in Eq.([31|) as 

^(P)oo ~ -Vaa(V = fyook+ip) (39) 

Substituting the latter representation into Eq.(|32j) and performing straightforward calcula- 
tions, we get that in the limit \i — > 00, the pairwise correlations in the asymptotic t — 00 
state obey 

Oaa(/i)oo ~ exp ( - — — (40) 
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where 

A - = ^kZ + kj (41) 

The approximation underlying the derivation of Eqs. (pfU|) and (|41|) is thus justified when 
or, in other words, when r chem = (AK+a^ + KJ) 1 > T diff = A 2 r /D. 
Hence, pairwise correlations in A particles distributions are long-ranged for microscopi- 
cally inhomogeneous reactions. A salient feature of the result in Eq. (|4()j) is that the charac- 
teristic decay length A corr of pairwise correlations depends on such "kinetic" parameter as 
D, the particles diffusion coefficient, which signifies that the asymptotic state achieved by re- 
versible diffusion-limited reactions with microscopically inhomogeneous elementary reaction 
rates at t = oo is not a true equilibrium state, but rather a non-equilibrium steady-state. 
We note finally that A corr is, of course, D-dependent for arbitrary relation between T c h em 
and Tdiff- Note also that for microscopically inhomogeneous reactions A corr does not vanish 
as D — > 0; in this case, the particles may perform random excursions in space just because 
in the break-up of B a pair of As is not forced to be born at their initial locations, at which 
they have entered the reaction forming the B particle. It is also interesting to remark that 
in the limit r c ^ em = [AK+cloo + KJ)~ l 3> r^// = A 2 r jD the correlation length A corr equals 
the distance travelled by a diffusive particle within its typical life-time between the reaction 
processes, i.e. the so-called Kuramoto length. Remarkably, precisely this len gth sets the 
scale of spatial correlations in out of equilibrium open chemical systems 

001 



B. "Contact reaction" approximation for microscopically inhomogeneous reactions. 

In this subsection we consider, within the framework of the decoupling procedure, the 
properties of the asymptotic t — > oo state for reversible diffusion-limited reactions in "contact 
approximation". This will allow us to illustrate the statement made in the Introduction 
that this approximation entails a somewhat ambiguous definition of the properties of the 
asymptotic t — oo state, as well as to present some explicit results on the corrections to the 
LMA in Eq.(j2J. 

Suppose that A;+(/x) and k-(fi) are some bell-shaped functions centered around their most 
probable values R and A, and the thicknesses of the distributions are negligibly small, such 
that the forward and the backward distance-dependent reaction constants can be deemed as 
delta-functions: k+(ji) = %\R)K + 5 d (\fi\- R) and = ■y~ 1 (X)K_5 d (\fi\ - A), where the 
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normalization factor 7d(i?) = A^R^ 1 , being the volume of a d-dimensional unit sphere, 
Ad = 2n d / 2 /T(d/2), and T(x) - the Gamma-function. 

For such a choice of the reaction rates, we get from Eq.(JB2} that 

<? aa {p)oo = F(p)(k^{p)b 00 - A;+(p)a^ - k + (p)a aa (\fi\ = R)co\ (42) 

which implies 

O"aa(/W)oo = J-{fJ>, A)&oo ~ J+(y, R) a lo ~ J +(^ R)o aa {\lA = R)oc, (43) 

where 

J-(H, A) = j^y d f dp e-^F(p)k4p), (44) 

and 

J+(/i ' R) = myJ dp e ~ lMF (p) k M- ( 45 ) 

Next, setting \p,\ = R in Eq. (|43jl . we obtain the following equation determining the value of 
the pairwise correlation function a aa (fj) at distance equal to the direct reaction radius: 

/. | m J-(|/i| = R, A)6oq - J + (|//| = R, g)c4 
^(W = i?)oo = i + J+ ( H=jR>jR ) ^ 

Consequently, taking advantage of Eqs. (|21|l and (|46|). we find that for reactions in contact 
approximation particles' mean concentrations &oo and obey: 

J = ^ (l + = ^, A) - f^-CH = A A)) . (47) 

Note now that setting i? = A, we get that </+(|/i| — R,R) — K + J(R) and = R, X = 

R) = K^J(R), which immediately implies that Eq. (|47|) simplifies to the conventional LMA in 
Eq. (J2J) and signifies that the chemical equilibrium state is a true thermodynamic equilibrium. 
If we, however, take into account dispersions around most probable values, which should be 
generally different for unimolecular backward and bimolecular forward elementary reactions, 
we would immediately obtain that the particles' mean concentrations do not obey the LMA 
and that the chemical equilibrium state is not a true equilibrium. The situation becomes 
even more ambiguous if the elementary reaction rates &+(//) and fc-(/i) are not simple bell- 
shaped, but more realistic complicated functions of /i having minima and maxima (see, 
e.g., Refs.0, Q)- As we have already remarked in the Introduction, for such k + (fi) and 
fc_(/i) the choice of R and A, at which the reactive boundary condition is imposed, has a 
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large degree of arbitrariness, but its consequences are grave: when one constrains himself 
to "contact approximation" and supposes, as it has been done in Refs.J^ and 2^, that 
two ill-defined reaction radii are equal to each other, one finds that there are no corrections 
to the LMA as t — > oo and that the CE is a true TES. On contrary, letting R and A be 
different on an arbitrarily small but fixed value, one gets that Eq. (|47|) is always different from 
Eq.(J2J), and that both and are dependent on such "kinetic" parameter as particles 
diffusion coefficient D. Which, in turn, implies that the asymptotic t = 00 state is not a 
true thermodynamic equilibrium, but rather a non-equilibrium steady-state. 

Consider next the explicit form of corrections to the LMA in Eq. (|47jl for microscopically 
inhomogeneous contact reactions in three dimensional systems. Here, we have for the Fourier 
transformed reaction constants: 



, . , sin(n_R) , , . . sin(pA) 

k + (p)=K + —^, and fc_(p)=ir_— ip, (48) 



and consequently, 



• 1 „ K+ . 1 „ , , K + f°° , / sm(pR) sin(pA) 

J + {\/i\ = R,R)-— ±J-(\n\=R,\) = —±- / pdpsm(pR)F(p) 1 



Vl ^' ' ' 2n 2 R J * * ^ ' ^' \ pR pX 

(49) 

Expanding next F(p) in Eq. ljHHj) into elementary fractions and performing the resulting 
integrals, we find that for sufficiently small R and A (A > R), 

J + (H =R,R)- §^-(M =R, A) « -£±-(1 - f ) (50) 

Substituting the last equation into Eq. (j4Tjl . we arrive eventually at the result presented 
in Eq.Q, which has been previously obtained in Ref. 21 1 in terms of a suitably extended 
heuristic Smoluchowski approach. In the next section we will show using a systematic 
diagrammatic expansion that, curiously enough, this equation is exact in the linear order in 
deviation from the equilibrium. 



IV. THE ASYMPTOTIC t = 00 STATE: EXACT RESULTS. 

In this section we will analyze behavior of the pairwise correlation function cr aa (p) enter- 
ing Eq. (J21|) in terms of a systematic diagrammatic expansion. In terms of this expansion, 
we will derive, in the linear with respect to the deviation from equilibrium approximation, 
an integral equation obeyed by cr aa (/V). We proceed to show then that our criterion of the 
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validity of the LMA and of the thermodynamic equilibrium, obtained in the previous section 
within the decoupling procedure, is exact. Moreover, we will define corrections to the LMA 
in Eq.pi for microscopically inhomogeneous reactions with contact and exponential reaction 
rates. 

We start by rewriting the Langevin equations (j!4j) in the following matricial form: 



where 5a = 5a(x,t) and 5(3 = 5/3(x,t) are local deviations from the stationary values, 
Eq. ffTBj) . C = ((x,t), e = e(p = O)^, the functionals are given by: 



while the linear operator M is defined in the Laplace- Fourier space of parameters s, p as: 



Note that in Eq.(|5^|). we have supposed, as in the previous section, that fc±(/i) depend only 
on the absolute value of /i and hence, that k±{p) = k±(—p). 

Denoting next the propagator G as an inverse of the matrix M, G = M , and taking 
the inverse Laplace transform, we find that the Fourier-transformed pair correlation function 
a aa (fi), Eq.(|THJ), obeys: 




(51) 




(52) 




(53) 



cr aa {p) = ([G aa {p,t)*{ Ul {p,t)-2e 5(p) + ({p,t) + 5a {p)5{t))] 2 ) 
+ ([GapM * {u 2 (p,t) + e 5(p) + 5f3 (p)5(t))] 2 ) 

+2 ([G aa (p, t) * ( Ul (p, t) - 2e 5{p) + C(p, t) + 5a {p)5{t))] 



(54) 



x [Ga/sM * (u 2 (p,t) + e 5(p) + 6/3o(p)6(t))]) 



where the symbol " * " denotes the time convolution, while 5ao and 5(3$ stand for the 
initial values of the Fourier spectra of fluctuations 5a (x, 0) and 5(3(x, 0). Note, however, 
that 5a and 5/3 give rise to exponentially decreasing terms and are insignificant for both 
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the asymptotic state and for the long-time kinetic behavior. Eq. (j54j) is complimented by 
the following equation determining the noise-noise correlation, Eq. (jl5|) . in time-momentum 
space: 

(C(p, t)C(~P, 0> = S(t - (e(p) - n(p)) (55) 
where e and f2 have been defined in the previous section. 

n n 

Now, following Refs. |27[ and [35| , we set up a systematic diagrammatic expansion of the 
pairwise correlation function <J aa (p), defined by Eqs.(p>fj) and (|53j). This expansion can be 
obtained either from the analysis of the action of the field theory describing the model 
or more directly by iterating Eq. (fo"^|) and collecting the different terms appearing during this 
procedure. The propagators and vertices corresponding to such an expansion are defined in 
figure H The two first kinds of vortices (of coordinance 3) arise because of the coupling of 
the quadratic term in ui and u<i either with 5a or with 5(3; the next one is due to the noise 
term ( and the last one to the e term. 

6a 6a 6a 6|3 6|3 6a 6|i 6|3 



time 



< 



5a 




5a 



5a 



5(3 




5a 5a _ 



5a 5a 




5a 



FIG. 1: Set of propagators (upper part) and vertices involved in the diagrammatic expansion of 
the correlation function a aa {jj). 




FIG. 2: Diagrammatic expansion of the correlation function a aa (p). 

A direct summation of these diagrams is impossible, of course. We will therefore deem 
this expansion as a perturbative development around the homogeneous situation, that is, an 
expansion around a configuration of the elementary reaction rates satisfying eQo)^ = 0. The 
small parameter in such an expansion will be the norm (supremum, for instance, denoted 
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as |e|) of the function e(p)oo. First order contribution can be obtained as follows: Note that 
a vortex e in figure [H corresponding to the term e in Eq. fjB^j) . is of order |e|; likewise, 
following Eq. (|55jl . one has that a vortex ( is also of order |e|. The first diagram of figure (J2J 
is therefore of order 1, the second of order 2, and all the following diagrams are of higher 
order in |e|. Consequently, in the linear order in deviation from the equilibrium, we obtain 
the following relation: 

-t 



dt'G» aa (p,t-ty e (p)-n(p) 



(56) 



where G aa {p, t) is the 0-order propagator, which can be readily obtained in an explicit form 
from Eq.(|53|): 

( K --<l~) r -(Dp 2 + q-)t + ( K --V + ) c -(Dy 2 + q + )t 



G° aa (p,t) = ^- —e ^ 1 * r + 

q + - q~ q~ — <T 



where 



K_ + 2 aoo (K + + fe + (p)) ± /2, 
K^ + 2a OD (K + + k + (p)')j 2 
- 8a 00 ( y K4K + + k + (p)) - 2k+(p)k.{p)y 

Turning next to the limit t —>■ oo, making use of Eq. (j21|) and denoting 

POO 

^°(p)= / dt'G° aa (P^ 
Jo 

we get the following system of implicit equations: 



,/\2 



(57) 



(58) 



(59) 



(60) 



(61) 



<? aa {p)co = F°{p) ^(p)oo - e(p) c 
e - ^(O)^ = 0, 

Note now that Eq. (|61|) has the same form as Eq. ()32|) . obtained in section 3 within the 
framework of the decoupling approximation. The only difference is in the form of the kernel 
J-°(p), which is a little bit more complicated than F(p) in Eq. (}3*2*|) but has exactly the same 
poles. Note, however, that this difference does not affect the criterium we have formulated in 
the previous section. One may immediately deduce from Eq.([61|) that cr a , Q ,(/x) 00 = 0, which 
implies that the LMA in Eq.Q is valid, if and only if e(p) equals 0, which holds only for 
systems with microscopically homogeneous reactions. This condition can be generalized to 
all orders and solves exactly our original problem in a very general framework for reactions 
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in Eq.Q: the chemical equilibrium is a true thermodynamic equilibrium and the LMA in 
Eq. (J2J) holds only for microscopically homogeneous diffusion-limited reversible reactions. If 
the reaction is not microscopically homogeneous, the asymptotic t = oo state (000,600) in 
the linear with respect to deviation from equilibrium approximation is described in a closed 
form by the system in Eqs. (j61j) . 



A. "Contact reactions" approximation revisited. 

In this subsection we turn back to the "contact reaction" case and re-examine it in 
within the framework of Eqs.([61jh which are exact in linear with respect to deviation from 
equilibrium approximation. As in section 3, we take the reaction rates in form of delta- 
functions: k+(p) = r y^ 1 (R)K + 6 d (\fi\ — R) and = 7 ( 7 1 (A)A'_<5 a! (|/i| — A) and suppose 
that A > R. 

^From Eq. (|57|) . we have that the Fourier-transformed 0-order propagator: 

G°-fe = , K K - „ e~ DpH + J K+a T K e- {Df? + 4A ' + °~ + K - )l (62) 

Integrating Eq.(|62|) over the time variable, we get that the kernel T°(p) in Eq. (J61j) is given 
explicitly by: 

1 ( Kl 16A>L 8K+K_a c 



(AK +aoo + K_) 2 \ 2D P 2 2(Dj? + 4A>oo + 2D P 2 + + K ~ 

(63) 

Note that this function of p, in the limit p — > 0, has exactly the same poles and shows 
exactly the same asymptotic behavior as F(p), Eq. (j33|) . obtained using an uncontrollable 
decoupling procedure. Solving next Eqs.lJoT]). we find that in three-dimensions particles 
mean concentrations a^ and 600 obey 



at K. 



00 



1 + s§s( 1 -f) +0 « A - fi ' 2 ) < 64 > 



K. 

This relation between the particles mean concentrations is obtained in an exact and con- 
trollable way. Quite surprisingly, in the linear with the respect to the difference (A — R) 
approximation, it coincide with Eq.fjHJ), obtained in Ref.j^jJ within the framework of a suit- 
ably extended heuristic Smoluchowski approach. 
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Similar to Eq.(|64|) results can be also obtained for low-dimensional systems. Performing 
straightforward calculations, we find then that and obey: 



K. 



1 + 



AuDR 



(X-R)+0[(X- R) 



and 



boo K. 

for two- and one-dimensional systems, respectively. 



l + ^±(\-R)+0[(X-R) 2 



(65) 



(66) 



B. Long-range, exponential reaction probabilities. 

Here we estimate corrections to the LMA for a particular example of long-range distance- 
dependent rates, characterized by an exponential dependence on the interparticle separation 
for the forward bimolecular reaction, and an exponential dependence on the radius of a 
geminate pair born in the elementary act of the unimolecular backward reaction. We focus 
on three-dimensional systems in which k + (fi) and fc_(A) obey: 

kM = iw^^ ^ M/i) = 8^ e ^ A (6?) 
For such a choice of reaction probabilities, we are not in position to solve our implicit 
equation (|61|) exactly and we have to resort to an approximate scheme. We thus develop 
an approximate approach, supposing that R and A are sufficiently small, such that the 
characteristic length A corr of correlations in particles distributions emerging due to ongoing 
microscopically inhomogeneous elementary reactions in Eq. (J67j) . is much greater than R and 
A. This is precisely the approximation used in section 3 to analyze behavior of correlations 
in general case of distance-dependent reactions. Under such an assumption, we may expect 
that o" aQ (/i)oo varies much slower than k + (/i) and fl(p)oo can be represented as in Eq. (j3T?j) . 
This implies that 

1 1 



a{p) 00 » b^K^ip) 

Evaluating next 



[l+p 2 X 2 ) 2 (l+p 2 R 2 ) 2 _ 



(68) 



n(ou w -K+ J tJ^Woo, (69) 
we get, using Eq.(|2*T|l. that particles mean concentrations ctco and 600 obey 
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which holds for R and A sufficiently small. Note that corrections to the LMA in case of 
long-range reactions rates in Eq. (j67|) appear to be smaller than in the "contact" case. 

Now, in order to verify that our approximation is self-consistent, we have to evaluate 
the form of pairwise correlation function cr aQ (/x) in /x-space. Performing the inverse Fourier 
transformation of (J aa {p)oc in Eq. (j68j) . we find that for sufficiently small R and A the pairwise 
correlations follow: 



AK + a^ + K_ 2 (4 K+aoo + K 



+C R e-^/ R + C x e~^/ X 

(71) 

where Cr and C\ are /i- independent constant and A corr has been defined in section 3, Eq.(JJTJ) 
in terms of the decoupling procedure. Therefore, our assumption underlying the derivation 
of the result in Eq.ffTTj) is justified when A corr 3> R, i.e., for sufficiently slow reactions such 
that the chemical time T c h em = (4-^+aoo + KJ)~ l is much larger than the time necessary 
to diffuse on distance R. Note also that the dependence of the correlation length on the 
diffusion coefficient is symptomatic of a non-equilibrium state. 



V. CONCLUSIONS. 



To conclude, in this paper we re-examined two fundamental concepts of classical chemical 
kinetics - the notion of "Chemical Equilibrium" and the "Law of Mass Action" - on example 
of diffusion-limited reversible A + A ^ B reactions. We have considered a rather general 
lattice model of such reactions, in which the elementary reaction rates are long-ranged 
and dependent on the instantaneous distance between particles. This model has been first 
analyzed by Zeldovich and Ovchinnikov [18] and pertains to chemical reactions between 
excited molecules or reactions involving transport of a proton or of an electron. It may be also 
viewed as a model of elementary reactions taking into account, albeit in an idealized fashion, 
different "microscopic" effects such as solvent structure, different angular orientations of 
reactive molecules, energy distributions and etc. In this model any pair of A particles, 
which perform standard random walks on sites of a d- dimensional hypercubic lattice may 
associate forming a B particle at any moment of time t with the rate k + (fi), where fi is the 
instantaneous distance separating these two particles. In the reverse reaction elementary 
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act, a mobile B particle may spontaneously dissociate with the rate fc-(A) giving birth to 
a geminate pair of As separated by a distance A. In our analysis we supposed that long- 
range rate k + (n) of bimolecular forward reaction and the rate k-(X), describing the birth of a 
geminate pair of A particles born at the distance A apart of each other within the elementary 
act of the unimolecular backward reaction, are arbitrary (integrable) functions of fi and A. 
n terms of a formally exact approach based on Gardiner's Poisson representation method 
jj, we have obtained exact non-linear Langevin equations describing the time evolution of 
complex-valued Poisson fields, whose mean values determine the A and B particles' mean 
concentrations. Solutions of these Langevin equations in the asymptotic t = oo state were 
obtained via a certain decoupling approximation and a systematic diagrammatic expansion. 

^From these solutions, which coincide in the leading order, we have deduced a general 
criterion determining the conditions when the classical LMA holds and when the asymptotic 
t = oo state is a true thermodynamic equilibrium. We have shown that this may only happen 
when the distance- dependent elementary reaction rates obey a very restrictive condition of 
microscopic homogeneity: the ratio /c + (/i)/fc_(/x) does not depend on \i for any /i. It seems 
that such a condition may be considered as apparently unrealistic since the bimolecular for- 
ward and unimolecular backward reactions are supported by different physical processes of 
classical and quantum origin. At present time, no general argument exists that it should be 
always the case. In case when is //-dependent, i.e. the elementary reactions 

are microscopically inhomogeneous, it appears that the detailed balance is broken, the LMA 
does not hold, particles' concentrations are spatially correlated and, remarkably, both cor- 
relation length and particles mean concentrations do depend on such kinetic parameter as 
the diffusion coefficient. This implies that the CE is not a true thermodynamic equilibrium 
but rather a non-equilibrium steady-state. Diffusion coefficient .D-dependent corrections to 
the LMA have been calculated explicitly in several particular cases. 

Consequently, for diffusion-limited reversible reactions the diffusional relaxation of the 
system is not fast enough to offset the perturbative effect of ongoing microscopically inhomo- 
geneous elementary reactions (breaking the detailed balance) even in the asymptotic t = oo 
state. This results in a non-vanishing current which modifies the structure of fluctuations 
and globally changes asymptotic concentrations. We emphasize that such a non-equilibrium 
steady-state emerges in a closed system with conserved overall concentration of particles 
and products, without any external inflow of particles. It might be also worthy to remark 
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that the diffusion-limited reactions provide thus a nice example of physical systems in which 
an arbitrarily small but finite difference (of classical or quantum origin) between the micro- 
scopic rates k+(ji) and entails a fundamental change in the asymptotic behavior of a 
many particle system. 

It is important to remark that our work does not contest any principle of thermodynam- 
ics, and in particular the zeroth principle. Indeed, in our model the system reaches a steady 
state, which can be described as the stationary point of a thermodynamical potential G, in 
agreement with the zeroth principle. Our work shows that G does depend on dynamical 
quantities such as the diffusion coefficient D, which implies that the entropic contribution 
to G can not be postulated a priori, but has to be defined using kinetic approaches taking 
into account dynamical effects. On the other hand, the LMA itself is not a law of ther- 
modynamics, but rather a law of an ideal gas, since it relies on the hypothesis of infinitely 
fast mixing of the reaction bath. Consequently, there is no surprise that the classic LMA 
does not hold for systems with microscopically inhomogeneous reactions and diffusion as the 
limiting transport process. 

Finally, we note that the analysis presented in our work may be extended in several direc- 
tions. In particular, the question of the corrections to the LMA may be addressed for other 
types of reaction schemes, including, e.g., reactions between excited molecules, in which case 
particles possess an intrinsic life-time. As well, one may expect that the effects observed here 
will become more pronounced for photochemical, ionization and electrochemical reactions 

H- 

The authors gratefully acknowledge helpful discussions with S.F.Burlatsky, S.Bratos, 
I.V.Gopich, M.Hankel, P.Hanggi, A.Lesne, K.Lindenberg, Z.Koza, M.Rausher, K.Seki, 
I.M.Sokolov and M.Tachiya. GO thanks the Alexander von Humboldt Foundation for the 
financial support via the Bessel Research Award. 



[1] J.Keizer, Statistical Thermodynamics of Nonequilibrium Processes, (Springer- Verlag, New 
York, 1987) 

[2] R.J.Silbey and R.A.Alberty, Physical Chemistry, (John Wiley and Sons, 2000) 



28 



[3] G.Wilemski and M.Fixman, J. Chem. Phys. 58, 4009 (1973); ibid. 60, 866, 878 (1974) 
[4] D.F.Calef and J.M.Deutch, Diffusion-Controlled Reactions, Ann. Rev. Phys. Chem. 34, 493 
(1983) 

[5] S.A.Rice, in: Diffusion-Limited Reactions, Comprehensive Chemical Kinetics Vol.25, eds.: 
C.H.Bamford, C.F.H.Tipper, and R.G.Compton (Elsevier, New York, 1985) 

[6] B.U.Felderhof and R.B.Jones, J. Chem. Phys. 103, 10201 (1995); J. Chem. Phys. 106, 954, 
967, 5006 (1997) 

[7] S.Lee and M.Karplus, J. Chem. Phys. 86, 1883 (1987) 

[8] K.L.Ivanov, N.N.Lukzen, A.B.Doktorov, and A.I.Burshtein, J. Chem. Phys. 114, 1754, 1763, 
5682 (2001) 

[9] C.W.Gardiner, Handbook of Stochastic Methods, (Springer, Berlin, 1990) 
[10] K.Tomita and H.Tomita, Prog. Theor. Phys. 51, 6 (1974) 
[11] Y. Kuramoto, Prog. Theor. Phys. 49, 1782 (1973); 52, 711 (1974) 
[12] G.Nicolis and M.Mansour, Phys. Re v. A 29, 2845 (1984) 

[13] K.Kitahara, K.Seki and S.Suzuki, J. Phys. Soc. Japan 59, 2309 (1990); K.Seki and K.Kitahara, 

J. Mol. Liq. 65/66, 293 (1995) 
[14] J.Gorecki, K.Kitahara, K.Yoshikawa, and I.Hanazaki, Physica A 211, 327 (1994) 
[15] J.Wakou and K.Kitahara, Physica A 281, 318 (2000) 

[16] F.C.Alcaraz, M.Droz, M.Hankel, and V.Rittenberg, Ann. Phys. 230, 250 (1994) 
[17] Ya.B.Zeldovich and A.A.Ovchinnikov, Sov. Phys. JETP Lett. 26, 440 (1977) 
[18] Ya.B.Zeldovich and A.A.Ovchinnikov, Sov. Phys. JETP 47, 829 (1978) 

[19] M.Tachiya, Proceedings of Ann. Meeting of Japanese Photochem. Assoc., Mie, Japan, October 
1980. 

[20] S.Redner and K.Kang, Phys. Rev. A 32, 435 (1985) 

[21] G.Oshanin, S.F.Burlatsky, and A.A.Ovchinnikov, Sov. J. Chem. Phys. 9, 588 (1991) 
[22] S.F.Burlatsky, A.A.Ovchinnikov, and G.Oshanin, Sov. Phys. JETP 68, 1153 (1989); J. Phys. 
A 22, 973, 977 (1989); Chem. Phys. 152, 13 (1991); G.Oshanin, Sov. J. Chem. Phys. 8, 395 
(1991); G.Oshanin, A.Mogutov, M.Moreau, and S.F.Burlatsky, J. Mol. Liquids 63, 175 (1995) 
[23] A.Szabo, J. Chem. Phys. 95, 2481 (1991); A.Szabo and R.Zwanzig, J. Stat. Phys. 65, 1057 
(1991) 

[24] J.Sung, K.J.Shin, and S.Lee, J. Chem. Phys. 107, 9418 (1997); M.Yang, S.Lee, and K.J.Shin, 



29 



J. Chem. Phys. 108, 9069 (1998) 
[25] W.Naumann, N.V.Shokhirev, and A.Szabo, Phys. Rev. Lett. 79, 3074 (1997) 
[26] I.V.Gopich, A.A.Ovchinnikov, and A.Szabo, Phys. Rev. Lett. 86, 922 (2001) 
[27] P.-A.Rey and J.Cardy, J. Phys. A 32, 1585 (1999) 

[28] I.V.Gopich and N.Agmon, Phys. Rev. Lett. 84, 2730 (2000); N.Agmon and I.V.Gopich, J. 

Chem. Phys. 112, 2863 (2000) 
[29] I.V.Gopich and A.Szabo, Chem. Phys. 284, 91 (2002) 

[30] D.Huppert, S.Y.Goldberg, A.Masad, and N.Agmon, Phys. Rev. Lett. 68, 3932 (1992); 

K.M.Solntsev, D. Huppert, and N.Agmon, J. Phys. Chem. A 105, 5868 (2001) 
[31] D.Pines and U. Pines, J. Chem. Phys. 115, 951 (2001) 

[32] N.Agmon and A.L.Edelstein, J. Chem. Phys. 100, 4181 (1994); A.L.Edelstein and N.Agmon, 

J. Mol. Liquids 64, 241 (1995); J. Phys. Chem. 99, 5389 (1995) 
[33] N.N.Bogolybov, Lectures on Quantum Statistic, (Gordon and Breach, New York, 1976) 



[34] S.Grossmann, J. C 
also Ch.6.5 in Ref. 



lem Phys. 65, 2007 (1976); J.Keizer, J. Chem. Phys. 67, 1473 (1977); see 
1 

[35] M.Droz and A.McKane, J. Phys. A 27, 467 (1994) 
[36] M.Doi, J. Phys. A 9, 1465, 1479 (1976) 
[37] L.Peliti, J. Physique 46, 1469 (1985) 

[38] S.H.Northrup and J.T.Hynes, Chem. Phys. Lett. 54, 244, 248 (1978); J. Chem. Phys. 69, 

5246, 5261 (1978); ibid. 71, 871, 884 (1979) 
[39] B. Gaveau, J. T. Hynes, R. Kapral, and M. Moreau, J. Stat. Phys. 56, 879, 894 (1989) 
[40] P.Hanggi, P.Talkner, and M.Borkovec, Reaction-Rate Theory: Fifty Years After Kramers, 

Rev. Mod. Phys. 62, 251 (1990) 
[41] W.Naumann and A.Molski, J. Chem. Phys. 103, 3474 (1995) 
[42] K.Seki, A.V.Barzykin, and M.Tachiya, J. Chem. Phys. 110, 7639 (1999) 
[43] G.J.Moro and M.G.Severin, J. Chem. Phys. 114, 4565 (2001) 

[44] S.Murata, S.Y.Matsuzaki, and M.Tachiya, J. Phys. Chem. 99, 5354 (1995); S.Murata and 

M.Tachiya, J. Chem. Phys. 100, 4064 (1996) 
[45] A.I.Burshtein, A.A.Zharikov, and N.V.Shokhirev, J. Chem. Phys. 96, 1951 (1992) 
[46] S.F.Burlatsky, Theor. Exp. Chem. 14, 483 (1978) 

[47] A.I.Burshtein and K.L.Ivanov, J. Phys. Chem. A 105, 3158 (2001); A.I.Burshtein, 



30 



A.A.Neufeld, and K.L.Ivanov, J. Chem. Phys. 115, 2652 (2001) 



